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BESCHRBIBUNG 

Verfidiren zum Segmenderen diner dreidimensionden Soi^^ 

Die Erfindung becrifit ein Verfehren ziim Segmencieren einer drcidimensionalen Struktur 
aus dnem dceidimensionaien Datensatz, der insbesondere tnedizinische Datenwerte 
enthalt, mit citiem deformierbaren Modell. Die Obetflachc des deformierbaren Modelk ist 
durch ein Netz mit Maschen gebildet, wobei die Masdhen Neczpunkte an der OberflXdie 
des Modelk miteinander verbinden. Die Erfindung bezicht sich aufierdem auf eine 
Bildverarbeimn^einriditung zur DurdiMirung dieses Ver^rens sowie auf dn 
Computerprogtamm zur Steuerung der Bildverarbeitungsdnricfatung. 

Ein Verfahren der eingangs genannten Art ist aus einem Aufsatz von Mclnem^j^ et al. 
„Deformabie modek in medical image analysis: asurvejr" in der Zeitschrift Medical Image 
Analysis, 1 (2):91-108, 1996, bekannt. Danadh wird ein deformierbares Moddl als dne 
elastisdie Oberflache dargestellt, die bspw, durch Minimierung einer gewichteten Summe 
aus einer incernen und einer externen sogenanncen Energie ihre Form und Position 
veiSndern kann. Die interne Energie kann dabei jede Grofie sein, die urn so kleiner wird, 
je weniger das Modell deformiert wird, wsthrend die externe Energie jede Grofie sein kann, 
die sich bei einer Annsiherung des defermierbaten Modelk an sogenannte 
I^dldatenpunktc, die wahrscbdjilidi auf der Obetflachc der zu segmenticrenden 
Struktur angeordnet sind, verkleinert. Ein Kandidatenpunkt kann bspw. aus einer 
voigegebenen Menge an Punkten im drddimensionalen Datensatz, deqenige Punkt sein, 
filr den die Wahrschdnlidikdt, auf der OberflSche der zu s^enrierenden Struktur 
an^ordnet zu sein, am grofiten ist. Ein Mafi fur diese Wahrschdnlidikdt ist die 
Strukturfunktion F^. Urn so grSjger der Wert der Strukturfunktion F(aO ist, desto grSfier 
ist die Wahrscheinlichkeit, dass der Punkt x im drddimensionalen Datensatz auf einfer 
Oberfiache der zu s^menderenden Struktur licgt. 

In bekannten Verfahren haben die Maschen z. B. die Form eines Dreiecks, dessen 
Edcpunkte durch drd benachbarte Netzpunkte an der Oberfiache des Moddk ddfiniett 
sind. Endang von Flachennormalen der Dreiecksflachen werden in dem dreidimensionalen 
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Dacensaiz Kaii<Mdatcnpunkte gcsucht, fUr die die StrukturfunJktion dn Maximum auiweist. 
Die Funktion FC«) karin dabei bspw. durch den Gtadienten der Datenwette an der Stdle x 
pto/izierc auf die FlSdbennomiale gegeben sdn. 

5 Nacfadem flir die Maschen des die Oberflache des Modells definierendcn Netzes 

Kandidatenpunkte etmittdt worden sind, erfolgt eine Neubctedinung det Netzpunkee des 
Modelis auf Basis der gefundenen Kandidatenpunkte. Dabd werden die Netzpunkte so 
berechnet, dafi fur die neu gefundenen Net^unkte die gewiditete Summe aus internet 
und externer Energie eiri Minimum eiiinimmt. Dariadi wird das Verfahren unter 
10 Zugrundelegung der neu beredineten Netzpunkee mehrfadi wiederholt, wobd sidi bd 
jeder Wiederholung das Moddl der Form der zu segmentierenden Struktur annahert. 
Dieses iterative Verfehren ist beendet. sobald dn Abbrudikriterium, z.B. eine besdmmte 
Anzahl an Wiederholungen, erfiillt ist. 

15 Naditd%istbddiesemVerfyiren,dassdncStrufcturiijnktionfbrdasgcsamte : 

defotmierbare Moded vcrwendet wird, so dass Strukturen, deren OberflSdie variierende 
Eigensdiaftcn aufweisen, nut unzureidjend s^endert werden kSnncn. Diese variierende 
Egensdiaft fcann bspw. ein Giadient der Datenweree sein, der filr untersdiiedKdie Teile 
der Oberfladie in en^engesetzte Riditungen zdgt. Bd andeten bekannten Vetfehrcn 

20 wird jeder Masche bzw. jeder Dreiedssfedie dne Struktutfunktion zugeordnet. was zu 
dnemhohenRcdienaufwandftihrt.AufierdemtretenbddiesemVeflfehrenoft 
Segmentienmgsfdiler auf. da bei dner Udnen Veisdiid>ung zwisdien dem defotmierbaten 
Moddl und der zu s^enticrcnden Struktur, nidit meht jede Strukturfunktion dem dazu 
geh6iigen Oberfladienbetddi der zu segmentierenden Struktur zugeordnet weidenkann. 

25 

Es ist Aufgabe der vorliegenden Erfmdung. dn Ver&hren anzugeben, bd dem diese 
Segmentierungsprobleme weniger ausgepragt sind. Diese Aufgabe wird gdSst durdi dn 
Vcifehren zum Segmentieren einer dreidimensionalen Struktur aus dnem 
drddimensionalen Daiensatz, der insbesondere medizinisdie Datenwerte endiSlt, mit den 
30 Sdiritten: 
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a) Bereitstellen eines dreidimensionalen, deformierbaren Modells, dessen Oberflacbe 
durch ein Netz mit Maschcn gebildet ist, die Necq>unkte an der Obcrftechc des 
Modells mitcinander vecbinden» wobei die IVIaschen in Gtuppai dngeteilt sind imd 
jeder Gtuppe eine Strukturfunkdon zugeoidnet ist, 

5 

b) Anordnen des Modells an der Stelle im dreidimensionalen Datensatz, an der sidh die 
zu sc^entierende Struktur befindec, 

c) Besdmmen je eines Kandidatenpunkces fur jede Masche mit Hilfe der 
10 Struktutfonkdon derjenigen Gruppe, deren Mi^lied die jeweilige Masche ist, 

d) Neuberedmen der Netzpunkte des ModeUs auf Basis der gefimdenen 
Kandidatenpunkte, 

15 e) Wiederholung der Sdiiitte c) und d) unter BerOcksichtigung der neu berechnecen 
Nec^unkte solange ein Abbnichl^iterium nicht erfMt ist. 

Im Ge^sisatz zu bekannten Vetfehren sind bei der Erfindung die Mascfien eines ModeUs 
in Gtuppen eingeteik und jeder Gmppe ist dne Strukturfunkdon zugeordnet; Dies 

20 erm6glichtdieBestimmungvonKandidatettpunkten,diemiteinerliohen 

Wahrscheinlichkdt auf dner Obetflachc der zu segmentierendcn Struktur lie^pn, auch bei 
einer Oberfladie, die stark in ihren Eigcnschaften -variiert, was zu einer verbessercen 
Segmenrierung fiihrt. Die variierende Eigenschaft fcann bspw., wie oben bereits erwahnt, 
die Richtung eines Gradicnten der Datenwerte sdn. Aufierdem erfordert dieses Ver&hten 

25 im Vergldch zur Zuordnung je dner Strukturfiuiktion zu jeder Masdie dnen geringeten 
. RechenaufWand, 

Bd den Ausgpstaltungpn gprnSg den AnsprUdien 2 und 3 wdsen die Maschen dner 
Gruppe ShnHdjeEigensdMiften auf. Da ahnliche Eigenschaften der Maschen auf ahnliche 
30 Eigpnschafien der entsprechcnden OberflSchenberdche der zu segmenderenden Struktur 
zuradaufiihren sind, ftihrt diese GtuppcnbUdung zu dner Differenzierung verschiedener 
Tdle der Struktur, z.B. der vetschiedenen Tdle dnes Organs oder eines Knochens. 
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Die Anspriiche 4 und 5 definieren bevorzugte Verfahren zur Besdmtnung von 
Kandidatenpimkten, 

5 Anspruch 6 defmiert eine Strukturfunktion, und Anspruch 7 beschreibc ein bevorzugtes 
Verfahren zur Neuberechnimg der Neczpunkte, das zu guten S^mentierungsergqbnissen 
fiihrc. 

Eine Bildvetarbeitun^einrichtung zur Durchfuhrimg des erfindungsgemafien Verlahrens 
10 ist in Anspruch 8 beschrieben. Anspruch 9 definiert ein Computeqprogtamm zur 
Steuerung der Bildverarbeicungseinrichtung nach Anspruch 8, 

Die Erfindung wird im Folgenden anhand der Zdchnungen naher erlautert. Es zeigen: 

15 Fig. 1 eine zur AusfUhrung des erfindungs^emSfien Verfahtens geeignece 
Biidverarbeitungsdnriditung in schemadscher Darstellung, 

Fig, 2 ein Ablau£liagcamm des erfindun^emafien Ver&hrens, 

20. Fig. 3 ein Ablaufdiagramm eines Vcrfahrens zur Gruppeneinteilung von Maschen eines 
defbrmierbaren ModeUs, 

Fig. 4 eine petspektivische Ansicht eines defbrmierbaren Modells einer linfcen 
Herzkammer, 

25 

Fig. 5 eine petspektivische Ansicht eines Femur-JVfodeUs und 

Fig. 6 zwei Helligkeitswerte-Vcrteilungen der mitderen Datenvektoren zweier Gruppen, 

30 Die in Fig. 1 dargescellte Bildverarbeitungseinrichtung umfesst einen Bildverarbeitungs- 
und Steuerprozessor 1 mit einem Speicher 2, in dem ein deformierbares Modell einer zu 
segmenrierenden Struktur und ein drcidimensionalcr Datensacz, der insbesondere 
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medizinische Datenwerce enchldc, gespeicherc sein koiuien. Der Bildvefarbeitungs- und. 
Sceuerprozessor 1 kann z.B. mittels eines Glad^rkabek 3 mit einer nicht nSher 
daigescellcen bildgebenden Einrichtung, wie einem Magnetresonanz- oder einem 
Computercomographiegerat, gekoppelt sein. Die von dem Bildverarbeittmgs- und 
5 Sceuerprozessor 1 s^mentierte Stniktut kann auf einem Monitor 4 wiedergegeben werden. 
Der Benutzer kann iiber eine Tascatur 5 oder durch andere, in Fig. 1 nicht naher 
dargesteilce Eingabeeinlieiten, auf den Bildverarbeicungs- und Sceuerprozessor 1 zugreifen. 

Fig. 2 zeigt den Ablauf eines Segmentierungsverfahrenst das mit der 
10 Bildverarbeitungseinrichtung nach Fig. 1 durchgefiihrc werden kann,: 

' Nach einer Iiutialisierung im Schritt 100 wird im Schricc 102 das deformierbare Modell 
einer zu segmentierenden Struktur, in dicsem Ausfiihrungsbeispiel ciner linken 
Hetzkammer, bereitgestellc. Die Oberflache des Modells ist dabei durch ein Netz mit 

15 Maschen gcbUdct, die Netzpunktc an der Oberflache des Modells miteinander verbinden. 
Die Maschen haben die Form von Dreiedcen, deren Eckpunkte durch Jeweik dtei 
benachbarte Netzpunkte an der Oberflache des Modells definiert sind. Anstelle eines in 
dicser Weise triangulierten Modelk k5nnte die Oberflache auch durch eine beUebige 
Polygonstrukcur definiert sein, bei dem die Maschen des Netzes nicht dreieckfermig sind, 

20 sondern eine von drei verscliiedene Anzahl an Netzpunkten unter Bildung dnes Polygons 
miteinander verbinden* 

Ein dcformierbares Modell kann aus mehreren TrainingsdatensStzen erzeugt wetden, in denen 
jeweils das gleiche Objekt durch einen Fachmann manuell s^mentiert wurde und die 

25 Oberflachen der segmentierten Objekte mit Netzen gleicher Topologie dargestellt werden. 
Gleiche Topologie bedeutet hier, dass jedes segmentierte Objekt die gleiche Anzahl an 
Maschen aufiveist und dass die gleiche Masche in uncerschiedlichen Trainingsdatensatzen die 
g^dchen Nachbarmaschen hat. Wenn also die Maschen in jedem Datensacz numerierc werden 
und bspw. im ctsten Datensatz die Masche Nr. 1 die Nachbarmaschen Nr. 2 bis Nr. 9 

30 aufwdst» so grenzt die Masche Nr. 1 auch in alien andcren Datensatzen an die Maschen Nr. 2 
bis Nr. 9. Eine Mittdung, im dnfechscen Fall eine arithmetische Mittdung* abet aUe 
segmentierten Objekte fUhrt zu einem Objekt. dass als deformierbares Modell bezeichnet 
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witd. 



Das erfindungsgemafie Vetfahren ist auch anwendbar, wenn das deformierbare ModeU in 
anderen Ausfiihrungsformen atif andere Art und Wdse etzeugt worden ist, z,B. durdi einfache 
5 Trianguladon eines Modells einer Enken Heizkammer. Ein hSufig angswandoes 
Triangularionsveriahren ist aus dem Ardkel ..Marching cubes: A high resolution 3d surfece 
construction algoridim.", W.E. Lorenwn und H.E. Cline, Computer Graphics, 21(3):163- 
169, 1987bckannt. 

10 Die Maschen des deformierbaren ModeUs werden nun im Schritt 104 in Gruppen eingeteilt, 
wobei das erfindungsgemafie Vetfehren mit jeder Methode zur Gruppeneinteilung 
ducdbfUhrbar isc. 

Eine Mdglicbkeit der GruppencinteUung ist in ciner cKten Ausaihrungsform in Fig. 3 
15 dargpsteflt Zunachst werden im Schritt 200 iVDatensatzc erzeugt. deten Datenwerte in 
diesem Ausfilhrungsbeispiel jeweils eine liake Hcrzkammer daistellen. Im Schritt 202 wird 
dann das deformierbare ModeB geti^ dem dngangs bcschtiebenen. bekannten Veifahren 
untcr Minimierung der gewichteten Summc der intemen und externea Energte an die 
Heizkammer in jedem Datensatz angcpasst. was zu A^JVTodeflen der Hetzkammer flihrt, die 
•20' alle die gleiche Topologie aufweisen. . - . 

ImSchritt204 wird eine Anzahl (?an Gruppen gpwahlt, indiedieDrdccksflScheneineeteUt 
werden. Wenn die Oberflache der zu segmenrierenden Struktur vetschiedeoe Bereicbe 
aufweist. die sich im dreidimensionalen Datensaiz edcennbar voneinander unteischeiden. so 

25 wird die Anzahl G bevorzugt so gewahlt. dass sie gleich der Anzahl an unteischiedlichen 
Oberilachenbereichen ist. Dieser erkennbare Unteischied kann bspw. ein Gradient der 
Datenwerte sain, der fUr einen Oberflachenbereich in das Innere der zu segmenrierenden 
Struktur gerichtet ist, und Bii einen anderen Oberflachenbereich nach aufien weist. Es ist 
bekannt, dass genau diese unterschiedlichen Gradientcnrichtungen bei der linken 

30 Heizkammer auftreten konnen. so dass in diesem Ausfiihr ungsbeispiel <?= 2 gewShlt wird. 

Im folgenden Schritt 206 wird jeder Dreiccksflache jedes Modells ein Merkmal zugeoidket. 
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Dieses Merkmal isc din Dacenvektor, der die Datenwerte (hier Helligkeicswecte) endiSlc, die 
auf der FIslchennormalen der Dreiecksflache liegen. Dabei verl^ufc die FlScbennoimale durdi 
den Miccelpunkt der Dreiecksflllche, und der Datenvdctor besteht bspw. aus. 10 Dacenwerten, 
die auf der Flachennormalen auf einer Seice der Dreiecks£&che liegen, aus einem Dacenwerc, 
5 der auf der DreiecksflSche liegt, und aus weiceren 10 Dacenwercen, die sich auf der aixderen 
Seite bzw. jenseits der Dreiecksilache befinden. In anderen AusfUhrungsfonnen konnten jeder 
Dreiecksflache auch andere Merkmale zugeordnet werden, z^B. der Gradient der Datenwerte 
an der Stelle des Flachentnitteipunktes projizicrt auf die Flachennormale. Auch mehr als ein 
Merkmal pro Dreiecksflache 5st niaglich. ' 

10 

Ein mittleres Merlcmal jeder Gruppe wird im Schritt 208 ermittelt. Dabei wird zunachst ein 
Merkmakraum, der alie Werte enthalc, die Merkmale, die den Dreiecksflachen zugeordnet 
wordensind, annehmen konnen» in C7Teilraume eingeteilt. Dann wird jeweils ein Mittelwert 
der Merkmale fiir jeden Teilraum bestimmt, und jeder Mitteiwert eines Teilraums wird als 
15 Gruppenmerkxnal definiert. In diesem Ausftihrungsbeispiel besteht der Merkmalstaum aus den 
Komponencen der Dacenvektoren, wobel jede Komponente Werte zwischen 0 und 511 
einnehmen kann. Der Merkmalstaum wird nun in zwd Tdbiaume au%eteilt. Die Merkmale 
des ersten Raumes sind Datenvektoren deren Komponenten Werte zwischen 0 und 255 
aufweisen, und die Merkmale des zweiten Raumes sind Datenvektox^n, deren Komponenten 
Werte zwischen 256 und 511 haben. Eine aridimetische Mittelung er^bt als mitderes 
Merkmal des ersten Teikaumes einen Datenvektor. dessen Komponenten den Wert 127.5 
haben. Der mittlere Datenvektor deszweiten Teilraumes weist entsprechend den Wert 383.5 
auf Somit ist das mitdere Merkmal der ersten Gruppe, also deren Gruppenmerkmal, der 
Datenvektor, dessen Komponenten den Wert 127.5 haben, und das entsprechende Merkmal 
der zweiten Gruppe ist der Datenvektor mit den Werten 383.5. 

Im Schritt 210 werden die zugeordneten Merkmale aller Dreiecksflachen mit Hilfe eines 
Ahnlichkeicsmafies mit alien Gruppenmerkmalen verglichen. Jede Dreiecksftelche wird 
dann der Gruppe zugeordnet. deren Gruppenmerkmal dem Merkmal der Dreiecksflache 
gemaK dem Ahnlichkeitsmafi am ahntichsten ist. Das Ahnlichkeitsmafi kann dabei jeder 
Ausdruck sein, der Merkmale der Dreiecksflachen bzw. der Gruppen miteinander 
vcrglcicht. Dies kann bspw. eine Korreladon der Datenvektoren zweier zu vei^eichender 
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Dreiecksftechen oder. wie in diescm Ausfuhrungsbeispiel, die Nonn der Difierenz zweier 
Datenvekcoren sein. 

Nachdem jede DreiecksflSche einer Gruppezugeordnet worden ist, wird im Schritt 212 ein 
5 neues mitderes Merkmal fiir jede Gruppe, also jeweik dn neues Gruppenmerkmal. ermittclt. 
Dazu wird uber die zugeordneten Mericmale aUer Dreiecke einer Gruppe gcmittelt. wobei in 

diesemAusfUhrungsbdspiduberdieDatenvekorenderDreieckearirfmiedsAgemit^^ 
Das Ergebnis dcr arithmetischen Mittelung ist das neue Gruppenmerkmal. 

10 I'nSchritt2l4wirduberpruft.obeinAbbruchkriteriumerfuUtist.DiesesAbbruchkrit^^ • 
fcann bspw. eine vorbesrimmte Anzahl an Iteradonen oder auch eine kaum vorhandene 
Anderung der Gruppenmerkmale sein. 1st das Abbrudikriterium etfullt. so witd mit Schritt 
21<5 fongcfahren. Ansonsten wird das Verfeliren mit dem Schritt 210 forigefiihrt. 

15 Da die ModeUe in den unterschiedlichen Datcnsatzen die gieiche Topologie aufWeisen. ist 
m jedem ModeU eine DreiecksflSche mit gldchen Nachbarschafisbeziehungen vorhanden. 
Dreiecke mit den gldchen Nachbarschafebeziehungen kSnnen nun in mdireren Gruppcn ■ 
Oder auch mdirmals m dner Gruppe endialten sdn. Diese Drdedce werden im Sdiritt 
216 der Gruppe zugeordnet in der diese Dreiecke am I^gsten enthalten sind. 
20 . . 

Aus den ATModdlen wird imSchritt218dndeformierbaresModdlduid.Mittdungeizeugt 
indem jeweik bspw. Qber die Posidonen aller Drdecksfl3chen, die die gfcichen 
Nachbarschafisbeziehungen aufWeisen, arithmerisch gemittdt wird. 
Ausfiihrlich ist dieses Verfahren in J'attem QassiBcarion and Scene Analysis". Duda and 
25 Hart, Wiley. 1973, besdirieben. wotauf hiermit vetwiesen wird. 

Die resulderende erste Gruppe 10. die zwcite Gruppe 12 und die Gruppengrenze 14 des ■ 
deformierbaren Moddls M der linken Herzkammer sind in Fig. 4 datgestdlt. Die Drdecke. 
die in diesem ModeU die erste Gruppe 10 bilden. entsprechen dem Berdch der linken 
30 Herzkammer. der die rechte Herzkammer berOhrt. D.h. durch die Gruppenemcdlung 
konnten verschiedene Teile des Herzens idendfizien werden. Diese Idenrifikadon einzdner 
Berddie der zu segmenderenden Struktur ist immer dann mSglich. wenn die im Datensatz • 
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erkennbaren Eigenschaften der Struktur, die die einzelnen Bereiche der Strufctur voneinaader 
unterscheiden, ak Merkmalezum Vergleich mittek eines AhnJichkeitsmafiesgenuczt werden. 
Die Idenrifizierung einzelaer Bereiche der zu segmenrierenden Stxuktur mittek der 
Gruppeneinteilung ermoglicht eitie Extraktion dieser Bereiche, so dass sie getreniit 
5 voneinander weiterverarbeitet werden kontien. 

Bin weiteres Beispiel dieser Idcntifizierung zeigt Fig* 5. Dort ist ein trianguliertes Femur- 
Modell F dargestellt, dessen Dreiecksflachen durch das oben besdhriebene. erfindung^emafie 
Verfehren in vier Gruppen 20, 22. 24, 26 aufgeteiit worden sind. Auch hier sind durch die 
10 Gruppeneinteilung Tcilc der zu segmendetenden Stniktur identiiizierc worden. So definicrt 
die Gruppe 22 den Femurkopf und die Gruppe 26 den gr6fieren Trochanter; 

Durch gedgnete Wahl der zu veigleichenden Merkmale und des AhnHdikeitsmdfies ist es abo 
m5glich, untetschiedlidhe Teile einer zu segmendetenden Struktur, wie bspw, eines Organs 
15 Oder eines Knochens, zu iden]d£izieren> die dann aus dem Datensatz extrahiert und individuell 
weiterverarbeitet werden konnen. 

In einer anderen AusRihrui^sform wird zur Einteilung der DreiecksiBachen des 
deformierbaren Modelk in Gruppen zunachst ein Drdeck willktirlich au^ewShlt, dessen 

20 Merkmal mit dem Merkmal eines angrenzenden, benachbarten Dreiedb veiglidien wird^ 
Merkmale sind hier die oben beschriebenen Datenvektoren. Das Verfkhrcn ist aberauch mit 
anderen Merkmalen ausfUhrbar. Verglichen werden die Dreiecksflachen mit dem ebenfalls 
bereits oben erlauterten Ahnlichkeitsmai?, wobei zwei Dreiecksflachen als ahnlich angpsehen 
werden, wenn das Mnlichkeitsmafi einen bestimmten Schwellwert iiber- oder linterschreitet. 

25 Dies kann bspw. ein Oberschreiten des Schwelhvertes 0,9 bei der Korrelation zweier 
Datenvektoren oder ein Untetschreiten des Schwellwertes 0,1 bei der Norm der DifFerenz 
zweier Datenvektoren sein. Sind die beiden Dreiecke ahnlich, so bilden die beiden Dreiecke 
eine erste Gruppe und ein oder mehrere Gruppenmerkmale werden durch Mittelung uber die 
einzelnen Merkmale der Dreiecke bestimmt- In diesem Ausflihrungsbeispiel kdnnten z.B. die 

30 einzelnen Komponcmen der Datenvektoren arithmetisdh gemitteit werden, so dass fiir diese 
Gruppe dn mitderer Datenvektor generiert wird. Dann wird genauso mit einem weiteren 
benachbarten Dreieck verfahren» wobd dann der Datenvektor dieses Dteiecks mit dem 
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mittleren Datenvektor der Gruppe verglichen wird. Dieses Verfahren wird fortgefohrt, bis 
keines der zu den Dreiecken der Gruppe benachbarten Dreiecke mehr dem mittleren 
Datenvektor der Gruppe ahnlich isc Dann wird ein weitetes, gruppenloses Dreieck 
ausgewahlt, mit dem unter Bildung einer zweiten Gruppe das eben beschriebene Verfahren 
wiederholt wird. Dieses Votgehensweise wird wiederholt> bis jedes Dteieck einer Gruppe 
zugeordnet worden isc. 

Eine weitere Ausfiihrungsform der Gruppeneinteilung ist die manuelle Einteilung durch einen 
' Experteii, der MascheA, die bezugiich eines oder mehrerer Merfcrnale der zusegmentierenden 
Struktur ahnlich sind, zusammenfasst, 

Im Schritt 106 wkd jeder Gruppe eine Strukturfunkdon zi^eordnet, wobei aJs 
Sttukturfiinktion jede Funkrion dicnen kann, die einen grofien Wert liefert» wenn ein 
Kandidatenpunkt mit einer grofien Wahrschcinlichkeit auf der OberflSche der zu 
segmentierenden Struktur UegL Da untetschiedliche Bereiche der Oberfilache der zu 
segmendetenden Struktur untetschiedliche Eigenschaften aufweisen konnen, ist cs sinnvoll, 
fiir diese Bereiche jeweils cine eigene Strukturfunkdon zu definieren. 

Ein Kandidatenpunkt kann eine Kandidaten-DreiecksflSche definieren. Hierbei ist die 
Kandidaten-DreiecksflSche eine zur Anpassung an eine zu segmenderende Struktur in 
Richtung der Flachennormalen einer Dreiecksflache des Modells versdiobene DteiecksfiSche. 
Der Kandidatenpunkt ist dann der Schnittpunkt der Kandidaten-Dteiecksflache mit der 
Flachennormalen. 

Zur Bestimmung der Strukturfunkrionen wird» wenn dies nicht bereits geschehen ist, uber die 
Merkmale der Dreiecksflachen einer Gruppe, also in diesem Ausfiihrungsbeispiel Qber die 
entsprechenden Datenvektoren, wie oben bereits beschrieben, gemittelt. Das resulrierende 
Gruppenmerkmal bzw, der mitderc Datenvektor definiert die Strukturfunktion der Gruppe. 
Den mitderen Datenvektor 30 der Gruppe 1 0 und den mitderen Datenvektor 32 der Gruppe 
12 zeigt Fig. 6. Dort ist auf der vertikalen Achse der Helligkeicswert und auf der horizontalen 
Achse ein Index, der die Datenwerte endang der Flachennormalen zahlt, au%etragen. 
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Die Strukturfunkdon ist in diesem Axisfiihrungsbeispiel ein Mafi jfiir die Cbeteinstimmung 
eines Merkmals eines Kandidatenpuaktes bzw. der dutch den Kandidatenpunkt definiertcn 
Kandidacen-Dreiecksflache mit dem entsprechenden Gruppenmerkmal, wobei als Mafi fur. die 
Obereinstimmungdas obenbeschriebeneAhnlichkeitsmafiverwendetwerden kann. Sokann 
5 die Strukturfunktion bspw. der Kehrwert der Norm der DifFerenz des Datenvektors einer 
Kandidaten-Dreiecksflache und des Datenvektors der entsprechenden Gruppe sein. 

Fttr jede Masche des defbrmierbaren Modells wird im Schri tt 108 eine Konfidenz mit einem 
Konfidenzmafi bestimmt. Die Konfidenz einer Masche gibt an, wie gut eine Masche in ihre 
10 jeweiltge Gruppe passt Bel einer hohen Konfidenz besteht eine grofie Sicherheit, dass diese 
Masche au%rund ihrer Merkmale wirklich der jewelligen Gruppe zugeordnet warden konnce. 
Bei einer gpringen Konfidenz ist die Zugehorigkeit zu einer Gruppe weniger sicher. 

Eine MogUchkeit, cin Konfidcnzmafi anzugeben, ist das oben beschriebenei^nlichkeit5ma& 
15 Sind die Merkmale eines .Dreiecks und einer Gruppe gemlUJ dem Ahnlichkeicsmafi Sehr 
ahnlicb, so ist auch die Konfidenz grofi. Das Resultat des AhnUchkeitsmafies konnte also 
direkt als Konfidenz genuczt warden. 

Ein weiteres Konfidenzmafi liefert die in Fig* 3 gezeigta Gruppaneintdlung. Dort wird im 
20 Schritt2l6furDreiecksflacJienmitgIeichenNadhbarsghafbbeEiehungOT^ 

diese Dreiecksflachen im Schritt 218 gemittelt worden ist> fiir jede Dreiecksflache, die 
Haufigkeit der Gruppenzugehori^eit bestimmt, Um so grcifJer diese Haufigkeit fiir eine 
Dreiecksflache isc, desto besser stimmen die Merkmale dieser Dreiecksflache mit den 
Gruppenmerkmalen uberein. Die Haufi^eit der Gruppenzugdiarigkeit kann also als 
25 Konfidenz verwendet werden. 

Das deformierbare Modell der linken Herzkammer wird im Folgenden zur automatischen 
Segmentierung einer linken Herzkammer in einem dieidimensionalen Datensatz gcnutzt. 

30 Dazu wird zun^chst im Schritt 1 10 das deformierbare Model! relativzur zu segmentierenden 
Struktur so ausgerichtet, dass die Oberflache des Modells moglichst gut mit den Kontuten der 
zu s^menrierenden Struktur ttbereinstimmc 
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Im Schritt 1 12 wird fiir jede Dreiecksflache der Obetflache des deformierbaren ModeUs jeweils 
ein Kandidatenpunkt x, ermittelc Der Kandidatenpunkt i, liegt dabei auf der 
Flachennormalen des Dreiecks tind auf einer Kandidaten-Dteiecksfteche. Die Position des 
5 Kattdidatcnpunktes x, auf der Flachennormalen wird nun so gcwahlt, dass die 
Struktutfiinkrion Fj(i,) der Gruppe. zu der die jeweilige Dreiecksflache gehorc, ein 

Maximum aufweist Die Indizes » undy bezeichnen hierbd verschiedene JGmdidatenpunkte 
brw. Gruppen, 

10 Nadidem fiir alle Drdedcsflachen der OberflScfae des deformierbaren Modells jeweils ein 
Kandidatenpunkt crmitteltworden ist, erfol^im Schritt 1 14 dneAnpassungder Obetflache 
des deformierbaren ModeUs an die gefimdenen Kandidatenpunkte. Dabd wild das Modellso 
deformiert, da& die sogenannte Enet^efunktton 



*5 +aE. 



20 



cxt 



(1) 



ein Minimum dnnimmt. Hierbd ist die externe Energie. die die DreiecksflSchen in 
Richtung der im Schritt 1 12 ermittdten Kandidatenpunkten bewegt, M^hrendE^ die interne 
Energie ist, die einer Deformation des Netzes entgcgenwirkt. Der Gewichtungsfektor 0 1st 
typischerweise so gew^hlt. dafi die bdden Enetgien in gidchcn Teilen zur Gesamteneigie 
beitragen. 

Als externe Energie kann jede Grofie herangezogen werden, die sich bd einer Annaherung des 
25 defonnierbarenModellsandiedurch dieBCandidatenpunktedefimerteStrukturverkleinert. In 
diesem Ausfiihrungsbdspid wird die externe Energie E„ durch folgende Gleichung 
besdurieben: 



30 



(2) 
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wobei uber samtliche Dreiecke summiert wird. Der Gewichtungsfakcor Wi berechnet rich nach ' 
folgender Gleichung: ^ .. . 

5 vf^,=max{0,F^(x,)}. 

(3) 

Gemafi Gleichung (3) werden Kandidatenpiuikce %\ deren Merkmale gut mit den 
Merkmalen der entsprcchenden Gruppe des Modelk ttbereinstimmen, starker gewichcet ak 
10 andere Kandidacenpunkte. 

In andeien Ausfuhrungsformen kdmiten die Gewichcun^fektoren auch die im Schritc 1 08 
ermittelte Konfidenz benicksichtigen« So kStinten Kandidatenpunkte, die einc geringere 
Konfidenz aufWeisen, weniger zur excemen Enet^e beitragen ab Kandidatenpunkte, fiir die 
15 eine groiSeres Konfidenz bestimtnt worden ist. Auch kdnnte die Sdiwankunsisbreite dec 
jeweiligen Gruppe berucksichtigt werden. Kandidatenpunkte dner Gruppe mit einer grofien 
Schwankungsbrcite, in denen also die Merkmale mehrerer Drciecksflachen von dem mitderen 
Merkmal der Gruppe stark abweichen, soUten mit cinem geringeren Gewicht berucksichtigt 
werden, ak Kandidatenpunkte einer Gruppe, die eine kleinere Schwankungsbreite aufweisc 

20 

Die interne Energie -£„,kann jede Grofie sein, die urn so kleiner wird, je weniger das Modell 
deformiert wird. In diesem Ausfiihrungsbeispiel berechnet sich interne Energie nach der 
Beziehung 

^ta. =S -X/)-jRAy)*. (4) 

25 Hierbei ist V die Anzahl der Preiecksfl2chen des Modells, s eln Skalierungsfaktor, R eine 
Rotationsmatrix und A^derAbstandzwischenzweiMittelpunkcenzweier Drciecksflachen des 
urspriinglichen, undeformierten Modelk Die Menge N(/) beinhaltet die Mitcelpunkte 
derjenigen Drciecksflachen, die zur Dreiecksflache mit dem Mitcelpunkt x,. benachbart sind. 
Bei der erstmaligen Minimierung der Energie gemafi Gleichung (1) ist a* = 1, und R ist die 
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Einheicsmacrix. 

Nachdem durch Minimierung der Energie gemaJU Gleichung (1) die Mittelpunkte Xj der 
Dreiecksflachen des deformierten Modells ermitcelt woxden sind, werden ein neuer 
5 Skaliemngsfaktor s und eine neue Rotatiommatrix R bestimmt, Dabei werden unter 
Berucksichtigung der aktuellen Mittelpunkte Xj der Skaliemngsfaktor und die Rotationsmacrix 
so gewaKIt, dass die interne Energie minimal wird In anderen Ausfiihruhgsformen konncen 
die iVIittelpunkte x^, der Skaiierungsfaktor s iind die Rotationsmatrix R gleichzeitig bestimmt 
werden, indem bei der Bestimmuhg des Mininiums der Gesamtehergie ^ aucK def 
10 Skaiierungsfaktor und die Rotationsmatrix variabel sind* 

Nachdem im Schritt 114 das deformierbare Modell an die zu segmentierende Struktur 
angepafic worden ist, wird im Schritt 1 16 geprUft, ob ein Abbruchkrirerium erfiillt ist. Dieses 
Abbruchkriterium feann z.B. eine vorbestimmte Anzahl an Wiederholungen der Schritte 1 12 
15 und 114^einodereinekaum nodi vorhandeneAnderuJtig der LagederNetzpunkte. Ist dieses 
Abbruchkriterium erfiillt, so endet das Segmentierungsver&hren im Schritt 118* Ist das 
Abbruchkriterium dagegen noch nicht erfUUt, so wird das Vetfahren mit dem Schritt 112 
fongefilhrt. 

20 Das nadi dem Segmentierungsverfehren erfaaltene Modell kann gespeichett und auf dem 
Momtor dargcstellt werden. Es stelit die Segmentierung der in dem dreidimensionalen 
Datensatz enthaltenen zu s^;mentierenden Struktur dar. 

Nachdem die Schritte 100 bis 108 einmal durchgefiihrt worden sind. ist dn deformierbaics 
25 Modell mic Maschen, die in Gruppen eingeteilt worden sind, und mit Strukturfunktionen, die 
den Gruppen zugeordnet worden sind, vorhanden. Daher kann bei der Segmentierung 
weiterer Strukturen auf dieses Modell zurQckgegriflfen werden, so dass die Schritte 104 bis 108 
entfallen konnen. 



30 
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BEZUGSZEICHENLISTE 

Veifahien zum Segmentieren eiaerdreidimensionalen Struktur 



M Deformierbates Modell.der link^n Herzkammer 

5 F Defonnierbares Femur-Modell 

1 Bildverarbeitungs- und Steuerprozessor 

2 Speicher 

3 GlasSiserkabel 

4 Monitor 
10 5 Tastatut 

10, 12 Gruppen von Dreiecksflachen des Modells der linken Herzkammer 

14 Grenziinie zwischen Gruppen 

20, 22, 24, 26 Gruppen von DrdecksfULclien des Femur-Modells. 

30, 32 Micderer Datenvektor einer Gruppe 
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PATENTANSPRtjrCHE 



1. Verfahten zum Segmenderea dner dreidimensioaalen Struktur aus einem. 
dreidimensjonalen Dacensaiz, der insbesondere medizinisdie Datmwerte enthalc, mit den 
Schritren: 



5 a) Bereitstellen eines dreidimensioaalen, defermlerbarcn Moddls, dessen OberBSxhe .. 

durdi dn Netz mit Ji^chen gebildet ist, die Netzpunkte an der OberflMdie des: 
Modells miteinander verbinden, wobd die Maschen in Gruppen eingeceiit sind und - • 
jeder Gruppe eine Scrukturfimkdon zi^eordnet ist, 

10 b) Anordnen des Modells an der Sceile im drddimenstonalen Datensatz, an der sidi die 
zu segmentierende Struktur befindet, 

c) Besdmmen je eines Kandidatenpunktes fiir jede Masche mit Hilfb der 
Strukcurfunktion derjenigen Gruppe, deren Mitglied die jeweilige Masdie ist, 

15 

d) Neuberedincn der Netzpunkte des Moddls auf Basis der gefimdenen 
Kandidatenpunkte, 

e) Wiederholung der Schritte c) und d) unter Berucksiditigung der ncu beredineten r 
20 Netzpunkte solange ein Abbruchkriterium nicht erfullt ist. 

2. Veriahren nadi Ansprudi 1, wobd der Schrin a) folgenden Schritte unifasst: 
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- Bereitstellen eines dreidimensionalen, deformierbaren Modells, dessen Oberflache 
durch ein Netz mit Maschen gebildet ist, die Netepunkte an der Oberflache des 
Modells miteinander verbinden, 

- Vergleichen jederMasche mit jeder angrenzenden, benachbarten Masche mit einem 
AhnlichkeitsmaB, 

- Zusammenfassen von gemSfl dem AhnlichkeitsmaB ahnUchen, benachbarten 
Masdien zu Gruppen, 

- Zuordnen je einer Strukturflmktion zu jeder Gruppe. 

3. Veifehren nach Aiwpruch 1, wobci der Schritt a) folgende Schritte umfasst: 

Bereitstellen von jeweils einem dreidimensionalen, deformierbaren Modell in 
mehreren DatensStzen, wobei alle Modelle die gleiche Anzahl an Maschen und die 
Maschen in jedem ModeU die gleichen Nachbaischaftsbeziehungen aufweisen, 

^ammenfassen der M asdien aller Modelle zu Gruppen, 

aiordnen der Maschaa mit gleichen Nadhtbarschaftsbeziehungen zu 
derjenigen 

Qmppe, in der diese Maschen am hftufigsten enthaltoi sind, 

- Zuordnenje einer Strukturfiinktionzu jeder Gruppe, 

- MittelungfiberdieModellederdrBidimensionalenDatensaizezurErmitthingeines 
deformierbaren Modells. 
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4. Verfahren nach Anspruch 1 , dadurcfa gdkennzeichnetv doss im Schrict c) aus einer Menge an 
Punkten derjetiige Punkt ak Kandidatenpunkt einer JMasche ennicceic wird> fUr den die 
Strukturfunktion ein Extremum, insbesondere ein Maximum aufweist. 

5 

. 5. Ver&hren nach Anspruch 1, dadurdi gekennzeidhmet, dass im Schrict c) ein 
Kandidatenpunkt einer Masche an der Stelle auf einer Normalen der dutch die Masche 
umsddossenen Fl^che ermittelt wird,, an der die Strukturfunktion ein Extremum, 
insbesondere ein Maximum aufweist. 

10 

6. Verfahren nach Anspruch 1, dadurch gekennzeichnec, dass der Wert der Strukturfimktion 

an einer besttmmten Stelle im dreidimensionalen Daten^atz ein MalS fUr die 
Wahrschdnlichkeit ist, dass diese bestimmte Stelle auf einer Oberfladie der zu 
segmentierenden Struktur angeordnet isc 

15 

7. Verfehren nach Anspruch 1, dadurch gekennzeichnet, dass im Schritt d) die Netzpunktfe 
unter Minimierung einer gewichtecen Summe aus extemer und interner Energie neu 
berechnet werden. 

20 8. Bildverarbeitungseinrichtung zur Durchfiihrung des Verfehrens nach Anspruch 1 mit 

einer Speichereinheit zum Speichem mindestens einos deformierbaren Modells, 
dessen Oberflache durch ein Netz mit Maschen gebildet ist, die Netzpunkte an der 
Oberflache des Modells miteinander verbinden, und zum Speichem mindestens 
25 eines dreidimensionalen Datensatzes, der insbesondere medizinische Datenwerte 

enfhat, 

einer Bildwiedergabeeinheit zur Wiedergabe einer zu segmentierenden bzw. 
segmentierten Struktur und des mindestens einen deformierbaren Modells, 

30 
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einerBeiechnungsemheitzurNeuborechnung derNetzpunkte desmindestens einen 
Modells, 

einer Anordmmgseinheit zum Anordnen des Modells, an der Stelle im 
dreidimensionalen Datensatz, an der sich die zu segmentierende Struktur befmdet, 

einer Steuereinheit zur Steuening der Speichereinheit, der Bildwiedergabeeinheit, 
der Berechnungseinheit und der Anordnungseinheit entsprechend den folgenden 
Scfaritten: 

a) Bereitstellen eines dreidimensionalen, deformierbaren 
Modells, dessen Obefflache durch ein Netz mit Maschen 
gebildet ist, die Netzpunkte an der OberflSehe des Modells 
miteinander verbinden, wobei die Maschen in Gruppen 
eingeteilt sind und jeder Gnippe eine Stnikturiiinktion 
aufweist, 

b) Anordnen des Modells an der Stelle im dreidimensionalen 
Datensatz, an der sich die zu segmentierende Struktur 
befindet, 

c) Bestimmen je eines Kandidatenpunktes fflr jede Masche 
mittels der Strukturfiraktion derjenigen Gnippe. deren 
Mitglied die jeweilige Masche ist, 

d) Neuberechnen der Netzpunkte des Modells auf Basis der 
gefundenen Kandidatenpunkte» 
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e) Wiederholimg <l«r Schritte c) und d) unter BerQcksichtigung der 
neu berechneten Netzpunfcte solange ein Abbruchfcriterium 
niditerftilltist 

9. Computerprogramm flir dne Steuereinheit zur Steuerung einer Speichereinheit, 
Bildwiedetgabednheit, einer Bercchnungseinhcit und einer Anoidnungseinheit 
Bildverarbeltungseinrichtimg zur Durchfiihrung des Vetfehrens nadi Anspruch 1 
entsprecfaend den folgenden Schrinen: 

a) BereitsteUen eines drddimensionalen. deformierbarcn Modells, dessen Obetflache 
durch ein Netz mit Maschen gebildet ist, die Netzpunkte an der Obeiflache des 
Modells mitdnander verbinden, wobd die Maschen in Gruppen dngpteUt sind und 
jede Gruppe eine Struknufunkdon auifweist. 

b) Anordnen des Modells an der Stdle im dieidimensionalen Datensatz, an der sich die 
zu s^menderende Strukmr befindet. 



c) Bestimmen je eines Kandidatenpunktcs ftir jcde Masche mitteb der Strukturfunkdi 
derjenigen Gruppe, deren Miiglied die jeweilige Masche ist. 



on 



d) Ncuberechnen der Netzpunkte des Moddk auf Basis der gefandeaen 
Kandidatetipunkte> 

e) Wiederholung det Schritte c) und d) unter BerUdcsichtigung der neu berechneten 
Netzpunkte solangp ein Abbruchkriterium nicht etfiiUt ist. 



PHDE030005 EP-P 



ZUSAMMENFASST TKfl 

Vetfahren zum Segmentteren einer dreidimensionalen Struktur 

Die Erfindung betrifft ein Verfehren zum Segmenrieren einer dreidimensionalen Struktur 
aus einem dreidimensionalen Datensatz, der insbesondere medizinische Datenwerte 
5 enthalt, mit Hilfe eines deformierbaren Modells. Die Oberflache des deformierbaren 
Modells ist durch ein Netz mit Maschen gebildet, wobei die Maschen Netzpunkte an der 
Oberflache des Modells miteinander verbinden und die Maschen in Gruppen eingeteilt 
sind. Jeder Gruppe ist eine Strukturfunkdon zugeordnet, die ein Mafi jflir die 
Wahrscheinlichkeit ist, daC ein Punkt des dreidimensionalen Datensatzes auf dec 

10 Oberflache der zu segmentierenden Struktur angeordnet ist. Nachdem das Modell an der 
Stelle im dreidimensionalen Datensatz angeordnet worden ist, an der sich die zu 
s^mcntierende Struktur befindet, wird fiir Jede Masche jeweils ein Kandidatenpunkt 
ermittelt, fttr den die Strukturfunktion der entsprechenden Gruppe ein Extremum 
aufweist. Die Netzpunkte des ModeUs wdcn auf Grundlage der Kandidatenpunkte neu 

15 berechnet 

Fig. 5 
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